The relationship of immune cells with autism spectrum disorder: a bidirectional Mendelian randomization study

Background Observational studies have indicated a correlation between immunological inflammation and the risk of autism spectrum disorder (ASD). However, the causal relationship between immunological inflammation and ASD remains uncertain. Methods Immunity-wide data sources were retrieved from the GWAS catalog. Genetic summary data on ASD were retrieved from two independent GWAS. We performed two independent bi-directional, two-sample Mendelian randomization (MR) analyses and a meta-analysis based on the two independent MR estimates to assess the causal relationship between ASD and immune cell signatures. Results We have discovered 26 potential correlations between genetic predisposition in the immunophenotypes and ASD. The meta-analysis of the two inverse variance weighted (IVW)-produced estimates provided further evidence supporting the potential causal relationship between immunophenotypes and ASD. Based on the findings of the reverse MR analysis, it was determined that there are two potential negative causal relationships between ASD and immunophenotypes. However, the meta-analysis of the two IVW-derived MR estimates indicated that immunophenotypes were not significantly influenced by ASD (OR = 0.87, 95% CI = 0.73 -1.03, P = 0.09; OR = 0.91, 95% CI = 0.81–1.01, P = 0.08). Conclusions This study expanded immune cell subtypes that were potentially causally associated with ASD risk as well as identified ASD-specific immune cell subtypes. The discovery has the potential to lead to earlier detection and more effective treatment techniques. Supplementary Information The online version contains supplementary material available at 10.1186/s12888-024-05927-5.

concern and a growing worldwide responsibility that has consequences for public health services [3].Multiple endeavors have been made to ascertain the cause of ASD, yet it still remains mostly unknown.The development of ASD is believed to be influenced by genetic, neurological, immunological, and environmental factors.Neuroimmunology is receiving increasing focus due to the potential connection between immune response dysregulation and abnormalities in neurodevelopment.Numerous studies have documented altered immune system function in people with ASD [4].
The long-held belief that the brain is an immune-privileged organ has shifted paradigms.In recent years, the interaction between the immune and nervous systems has been well described in a variety of pathological cases, ranging from autoimmune diseases like multiple sclerosis to more traditional cognitive dysfunctions like Alzheimer's disease [5].Numerous studies have shown that people with ASD have a different immunological response.Evidence of changes in the functioning of the central and peripheral immune systems suggests that a subgroup of people with ASD exhibit immunological dysregulation [6].Furthermore, children with ASD have a higher frequency of immune-related comorbidities such as autoimmune disorders, allergies, and psoriasis when compared to healthy controls [7].Immune system changes include inappropriate immune cell activation, the production of autoantibodies, cytokine and chemokine imbalances, and increased permeability of the blood-brain barrier [8].Individuals with ASD often exhibit lymphocyte immunomodulation (T lymphocytes, B lymphocytes, natural killer (NK) cells, and dendritic cells).The onset of maternal autoimmune illnesses during pregnancy might significantly influence the likelihood of ASD in the child as well [9].On the other hand, the current research has shown contradictory conclusions regarding the connection between immunological inflammation and ASD [10].This might be attributed to factors such as small sample numbers, errors in study design, and the presence of confounding variables that were not well addressed in these investigations.ASD has a high level of complexity and heterogeneity.The topic of whether immunological dysregulation is a fundamental etiology or a later outcome remains unresolved [8].
Mendelian randomization (MR) is a method that mimics a randomized controlled trial and utilizes genetic information as an instrumental variable to assess exposure.Because allele allocation occurs during meiosis and is not influenced by environmental exposures, MR analysis is considered to be less affected by underlying confounding factors and may be used to infer causal relationships, provided the statistical assumptions are satisfied [11].However, to our knowledge, no MR study has been conducted to evaluate the correlation between immune cells and ASD.To further investigate the causal roles played by peripheral immunity in ASD risk, the present research used two-sample Mendelian randomization (MR) studies to evaluate the possible causal association between immune cells and ASD.

Study design
This research made use of publicly accessible data.All original investigations received ethical clearance and informed permission from their subjects.To investigate the causal connection between ASD and 731 immune cell signatures (7 groups), we conducted two separate bidirectional, two-sample MR studies and a meta-analysis based on the two independent MR estimations.The MR technique involves three steps: (i) identification of instrumental variables (IVs) via the use of genetic variations; (ii) evaluation of causal estimations; and (iii) conducting sensitivity studies.Figure 1 displays the conceptual representations of our investigation.

Immunity-wide GWAS data sources
The GWAS summary data for each immunological feature may be accessed publicly via the GWAS Catalog.The accession codes range from GCST90001391 to GCST90002121 (https:// www.ebi.ac.uk/ gwas/ studi es) [12].A comprehensive set of 731 immunophenotypes was analyzed, which included several measurements such as absolute cell counts (n = 118), median fluorescence intensities (MFI) indicating surface antigen levels (n = 389), morphological parameters (MP) (n = 32), and relative cell counts (n = 192).More precisely, the MFI, AC, and RC characteristics include B cells, CDCs, fully developed T cells, monocytes, myeloid cells, TBNK (T cells, B cells, natural killer cells), and Treg panels.On the other hand, the MP characteristic includes CDC and TBNK panels.The first genome-wide association study (GWAS) on immunological variables included data from 3,757 people of European descent, with no cohorts overlapping.Around 22 million single nucleotide polymorphisms (SNPs) were genotyped using high-density arrays.These SNPs were imputed using a Sardinian sequencebased reference panel [13].Associations were then assessed, taking into account variables such as sex, age, and age squared (Table 1).

GWAS data sources for ASD
We discovered a grand total of 7 ASD datasets published between the years 2013 and 2021 while doing a search for ASD datasets in the IEU OPEN GWAS PROJECT database.In order to prevent the issue of having duplicate inclusion samples, we opted for two datasets, namely ieu-a-1184 and ieu-a-806, which were both released in the same year.Simultaneously, we considered factors such as race, sample size, and the quantity of SNPs.The datasets were both published in 2015 and consist of individuals of European origin.The first dataset has a sample size of 10,610, while the second dataset has a sample size of 10,263.The SNPs identified were 9,499,589 and 9,499,590, respectively (Table 1).

Selection of instrumental tools
A series of quality control measures were conducted to choose appropriate genetic instrumental instruments [14].Based on current research [15,16], we examined the single nucleotide polymorphisms (SNPs) that are linked to each immune trait at a significance level of p < 5 × 10 -5 , with a clumping window size greater than 10,000 kb and a linkage disequilibrium level of r 2 < 0.001.In the case of ASD, we modified the significance threshold to 5 × 10 -8 .The F statistic was used to verify the robust correlation between IVs and exposure.A value of the F statistic over 10 was often deemed to satisfy the criteria for a strong connection [17].

Statistical analysis
In this work, we conducted two separate bidirectional, two-sample MR studies.The research used four MR approaches to investigate the causal relationships between 731 immunophenotypes and ASD: inverse variance weighted (IVW), MR-Egger method, weighted model, and weighted median method.The IVW methodology was chosen as the primary method because it used a weighted linear regression model.The Cochran's  Q statistic and its accompanying P-values were used to assess the heterogeneity among the chosen IVs.In the event that the null hypothesis is rejected, the fixed-effects IVW method is substituted by the random-effects IVW method [18].We used the MR-Egger regression intercept test to detect the presence of horizontal pleiotropy.
In addition, the MR pleiotropy residual sum approach was used to eliminate any horizontal pleiotropic outliers that might significantly impact the estimate findings.Furthermore, scatter plots and funnel plots were used.
Ultimately, we used a fixed-effect model to conduct a meta-analysis on the two MR estimations produced from the IVW-MRE method.The statistical analyses were conducted using the R software (version 4.0.3).The MR analysis was conducted using the TwoSampleMR programs [19].The TwoSampleMR package enables the execution of Mendelian randomization using GWAS summary data.It automatically retrieves data from the IEU GWAS database and offers a wide range of methods for conducting the analysis.In order to tackle the problem of multiple testing, we used a Bonferroni-corrected significance criterion [20,21].This threshold was determined by dividing 0.05 by the total number of immunophenotypes (731), resulting in a value of 0.000068.P values ranging from 0.000068 to 0.05 were deemed to suggest potential causal relationships between the exposures and the outcomes.Causal connections between the exposures and the outcomes were judged significant if the P-value was less than 0.000068.

Investigation into the causative impact of immunophenotypes on ASD
IVW was selected as the principal MR method.According to the analysis of the ieu-a-1184 dataset, there are sixteen immunophenotypes that have a beneficial effect on ASD, whereas there are eighteen immunophenotypes that are associated with an elevated susceptibility to ASD (Fig. 2A).The results of the Bonferroni-adjusted analysis show a strong causal relationship between ASD and CD3-lymphocyte %leukocyte (P = 0.00002), additionally, there is a potential causal relationship between ASD and the other 33 immune phenotypes (Table 2).The analysis of the ieu-a-806 dataset revealed that fifteen immunophenotypes have a protective impact on ASD, whereas sixteen immunophenotypes are linked to an increased vulnerability to ASD (Fig. 2B).The Bonferroniadjusted analysis indicates a substantial causative connection between ASD and CD3-lymphocyte%leukocyte (P = 0.00001).ASD may have a potential causal  3).

Investigation into the causative impact of ASD development on immunophenotypes
In order to investigate the direct impacts of ASD on immunophenotypes, a two-sample MR study was conducted, with the IVW method used as the primary analytical approach.Through the GWAS of ieu-a-1184, the results of the Bonferroni-adjusted analysis revealed 13 immunophenotypes with potential indications (Fig. 5A and Table 4).GWAS analysis of ieu-a-806 identified 21 immunophenotypes with potential correlations based on Bonferroni-adjusted analysis(Fig.5B and Table 5).By identifying the common features between the two groups, we discovered a potential causal relationship between a greater risk of ASD and a lower level of circulating CD19 on IgD + CD38-unswitched memory B cells and CM CD4 + %CD4 + T cells.The meta-analysis of the two IVW-derived MR estimates resulted in an o OR of 0.87 (95% [CI] 0.73-1.03,P = 0.09) for CD19 on IgD + CD38-unswitched memory B cells and 0.91 (95% [CI] 0.81-1.01,P = 0.08) for CM CD4 + %CD4 + T cells (Fig. 6).

Discussion
Numerous efforts have been made to establish the etiology of ASD; nonetheless, it remains mainly mysterious.The development of ASD is believed to be influenced by genetic, neurological, immunological, and environmental factors.Neuroimmunology is receiving increasing focus due to the potential connection between immune response dysregulation and abnormalities in neurodevelopment.Multiple studies have shown altered immune system function in patients with ASD [22].In this study, we used a bi-directional, two-sample Mendelian randomization method to establish that 26 specific immunophenotypes had potential causal effects on ASD across four distinct immunological characteristics (MFI, RC, AC, and MP).This work is the first effort to uncover the cause-and-effect connection between 731 immune cell characteristics and ASD using a genetic method that relies on GWA summary data.T-lymphocytes have a crucial function in facilitating adaptive immunity through the secretion of cytokines, which help regulate immune responses at the local level    [27].Furthermore, a considerable proportion of individuals with ASD have impairments in CD4 + and CD25 + Tregs, resulting in the activation of the immune system against self-tissues in a specific subset of these patients [28].In addition, a recent meta-analysis has shown significant abnormalities in CD4 + lymphocytes, namely a decrease in Tregs and an increase in Th17 cells, in individuals with ASD [27].Treg cells have a crucial function in maintaining self-tolerance by restraining autoimmune reactions via the suppression of proinflammatory mechanisms [29].Several investigations have shown abnormalities in lymphocytes and imbalances among different types of lymphocyte subgroups.The most consistent results in individuals with ASD include reduced responsiveness to stimulation [30,31], aberrant activation [32], an  imbalanced ratio of T helper and suppressor cells [33], decreased activity of Th cells [34], and a systemic deficiency of Tregs [35].The findings as a whole indicate that a significant percentage of patients with ASD have changes in the function of their lymphocytes, particularly in T-cell subpopulations.These findings, in conjunction with cytokine irregularities, provide a more comprehensive perspective on a potential underlying cause for the reported anomalies in ASD.Natural killer (NK) cells make up around 15% of the lymphocytes found in the bloodstream and have a crucial function in the innate immune system [36].These cells are distinguished by the absence of CD3 surface antigen and the presence of CD56, and their activity is carried out by the secretion of immunomodulatory cytokines like IFN-γ, tumor necrosis factor-alpha (TNF-α), and IL-10.In addition, they possess cytolytic activity and facilitate cellular cytotoxicity and immunological surveillance by interacting with dendritic cells [37,38].Discrepancies in the activation and inhibition states of some factors may contribute to the development of autoimmune disorders.However, the precise processes responsible for this phenomenon have not been completely elucidated [39].NK cells have been implicated in the development of neurological illnesses such as multiple sclerosis [40], schizophrenia [41], Tourette syndrome [42], and Rett syndrome [43].The impact of NK cells on ADS is a subject of debate in research.Two prior studies [44,45] found that autistic people had larger absolute numbers of NK cells in their peripheral blood.While another study found that a large number of people with ASD had a decrease in their NK cell count [46], our study discovered that increasing the fraction of CD3-lymphocytes reduced the incidence of ASD.The selection and quantification of NK cells were based only on the use of CD56 and CD3 markers.CD56 and CD57 are markers that can be utilized to categorize many discrete subsets of NK cells.The transcriptional signature, phenotypic characteristics, and functional abilities of human CD57 + NK cells differ from those of CD57 + -NK cells [47,48].Siniscalco et al. discovered a decline in CD57 + CD3-numbers among individuals with ASD, although D56 + CD3-counts remained unchanged but within the normal range [46].This suggests that a certain group of non-T lymphocytes with NK activity may play a crucial role in ASD.
A correlation between certain alleles of human leukocyte antigens (HLA) and autoimmune disorders has been discovered.Multiple studies have established connections between HLA and ASD, demonstrating that autistic children have a greater occurrence of the HLA-DRB1*11 allele and a reduced occurrence of the HLA-DRB1*03 allele [49].Multiple researchers have documented a correlation between HLA and ASD in various populations, such as Caucasians [50] and Chinese [51].Guerini et al. discovered an intriguing correlation between HLA-G polymorphism and ASD, which is likely due to prenatal immunological activation [52].The presence of HLA-DR4 in mothers has been identified as a risk factor for ASD in their children [53].When investigating the immunopathology of the disease, it is recommended to take into account genetic variations in the HLA region [54].Our MR analysis indicated a potential causal association between HLA DR expression on plasmacytoid dendritic cells, the percentage of monocytes expressing CD62L-HLA DR + + , HLA DR expression on CD33-HLA DR + cells, and CD45 expression on HLA DR + CD4 + TBNK cells, and ASD.This demonstrates the importance of HLA in the pathophysiology of ASD.
In individuals with ASD, there is an occurrence of abnormal numbers or activation of certain innate immune cells.This includes a notable increase in the quantity of peripheral myeloid dendritic cells (DCs) [55].The alterations in the DC population might be attributed to an augmented process of monocyte differentiation into dendritic cells.Additional research has shown that peripheral myeloid DCs obtained from children with ASD have a higher level of surface CD80/CD86 co-stimulatory molecules compared to typically developing children [56].DCs variations were shown to be associated with differences in amygdala volume, gastrointestinal problems, and impaired behaviors [55].Our study found that the presence of CD80 on myeloid DC and HLA DR on plasmacytoid DC is associated with an increased risk of ASD.Conversely, the presence of CD62L on CD62L + DC, CD80 on plasmacytoid DC, and CD80 on CD62L + plasmacytoid DC is associated with a decreased risk of ASD.This indicates that different subtypes of dendritic cells have differing effects on the risk of ASD.
Heuer et al. examined the roles of B cells and noted a reduction in the overall levels of both IgM and IgG in the peripheral blood of children with ASD, as compared to typically developing individuals [57].They subsequently examine whether the decreased amounts of IgG and IgM in the plasma were caused by impaired formation, activation, or activity of B cells.Heuer et al. demonstrated that there were no discernible disparities in the quantity of B memory cells.These findings suggest that the reduction of immunoglobulins (IgG) in individuals with ASD is not caused by B cell malfunction but rather is influenced by the participation of many immune cells [58].Furthermore, based on a meta-analysis comprising five investigations, no statistically significant variations in B cell counts were seen between persons with ASD and normally developing controls [59].Our findings indicate the presence of eight distinct B cell subtypes, with five showing a positive causal relationship and three showing a reverse causal association with ASD.Additional population-based and experimental studies are necessary to fully elucidate the intricate connection between specific subtypes of B lymphocytes and ASD in the future.
Naturally, this study has certain constraints.Initially, there is a possibility of selection bias due to the fact that the immune cell samples were mostly from individuals of Sardinian descent, whereas the samples of individuals with ASD were exclusively from European ethnicities.The Sardinian population, used for the GWAS of immune cell features, is genetically distinct as a result of its geographical isolation from mainland Europeans.Orru et al. (2020) discovered that out of the 122 genetic signals related to immune cell features, 16 (13%) had low frequencies among Europeans from the 1000 Genome Project, but high rates among Sardinians [12].Furthermore, ASD is predominantly diagnosed in males, with a male-to-female ratio of up to 4:1.Our study did not distinguish between genders.In order to accurately interpret our findings, it is imperative to assess potential disparities between the outcomes when applied only to male or female populations.
In conclusion, we thoroughly investigated the causal relationship between the 731 immunophenotype and ASD.Our results suggest that there are eleven potential positive causal directions and fifteen potential reverse causal directions for ASD.The study indicates that specific types of immune cells and genetic predispositions could serve as biomarkers for the likelihood of developing ASD.This discovery has the potential to enable earlier detection and more efficient treatment strategies.

Fig. 1
Fig. 1 The diagrammatic depictions of our investigation.The IVW model refers to the inverse variance weighted model

Fig. 2
Fig. 2 Forest plots illustrate the causal associations between immune cell characteristics and ASD by using IVW methods.A Estimates based on the ieu-a-1184 dataset; B Estimates based on the ieu-a-806 dataset

Fig. 3 Fig. 4
Fig.3The meta-analysis of the two IVW-derived MR estimates revealed fifteen positive causal directions associated with ASD.The MR results are based on the GWAS of the ieu-a-1184 dataset and the ieu-a-806 dataset

Fig. 5
Fig.5shows the causal associations between ASD and immune cell characteristics by using the IVW method.A Estimates based on the ieu-a-805 dataset; B Estimates based on the ieu-a-1148 dataset

[ 23 ]
. T-lymphocytes are comprised of two subsets: CD8 cytotoxic T cells and CD4 helper T cells.CD4 T cells that have been activated undergo cellular division and transform into distinct subsets of T-helper (Th) cells.The subtypes mentioned are Th1, Th2, Th17, and regulatory T (Treg) cells[24].Treg cells have a role in avoiding the occurrence of autoimmune diseases and maintaining the body's ability to tolerate itself by releasing the antiinflammatory cytokines IL-10 and TGF-β[25].IL-10 is produced by Th2 and Treg cells, along with various other immune cells, in order to restrict the excessive release of pro-inflammatory cytokines and chemokines during infection and reduce tissue damage[26].Our study discovered that there is a negative correlation between the risk of ASD and the proportion of CD3 on CD4 + regulatory T cells decreases.An analysis of immune cells in the blood of individuals with ASD has shown a significant decrease in Treg compared to healthy individuals, as demonstrated by De Giacomo et al. in 2021

Fig. 6
Fig. 6 Meta-analysis estimates of MR results from the IVW model based on the GWAS of the